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Detecting the Cosmic Dipole Anisotropy in Large-Scale Radio Surveys 

Pronefield Crawford^ 

ABSTRACT 

The detection of a dipole anisotropy in the sky distribution of sources in large-scale 
radio surveys can be used to constrain the magnitude and direction of our local mo- 
tion with respect to an isotropically distributed extragalactic radio source population. 
Such a population is predicted to be present at cosmological redshifts in an isotropically 
expanding universe. The extragalactic radio source population is observed to have a 
median redshift of z ~ 1, a much later epoch than the cosmic microwave background 
(z ~ 1100). I consider the detectability of a velocity dipole anisotropy in radio sur- 
veys having a finite number of source counts. The statistical significance of a velocity 
dipole detection from radio source counts is also discussed in detail. I find that existing 
large-scale radio survey catalogs do not have a sufficient number of sources to detect 
the expected velocity dipole with statistical significance, even if survey masking and 
flux calibration complications can be completely eliminated (i.e., if both the surveys 
and observing instruments are perfect). However, a dipole anisotropy should be easily 
detectable in future radio surveys planned with next-generation radio facilities, such 
as the Low Frequency Array and the Square Kilometer Array; tight constraints on the 
dipole magnitude and direction should be possible if flux calibration problems can be 
sufficiently minimized or corrected and contamination from local sources eliminated. 

Subject headings: methods: analytical - catalogs - large-scale structure of universe 

1. Introduction and Motivation 

The measurement of a dipole anisotropy in large-scale surveys can be used to probe the dis- 
tribution of matter at different distances and constrain our local motion with respect to large-scale 
mass distributions in the universe. There have been a number of attempts to measure the dipole 
anisotropy in different surveys. The dipole anisotropy in the cosmic microwave background at red- 



shift z ~ 1100, probed by the Cosmic Backgro und Explorer (ISmoot et aLlllOQ^ l. and more recently 



by the Wilkinson Microwave Anisotropy Probe (jBennett et al.ll2003l ). was clearly detected, and tight 
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constraints on our local motion with respect to this background were m ade with these mea sure- 
ments. Analysis of the X-ray background at intermediate redshifts {z < 4) (jBoughn et al.ll2002l ) and 
the matter distribution at local distances (the Local Group), probed at in frared wavelengths using 
data from the Infrared Astronomical Satellite (IRAS) (IWebster et al.lll997l) . have resulted in dipole 
detections as well. A marginal detection wa s also reported bvlBlake &: Walll (j2002l ) using radio data 
from the NRAO VLA Sky Survey (NVSS) (jCondon et alJ ll998l). In Section 4 2, I discuss this radio 
survey specifically in the context of a dipole measurement. iBaleisis et al.l (119981) refer to a large 



number of previo u s dipole investigations using infrared galaxy catal ogs (e.g.. lMeiksin &: Davis 



Yahil et al. 19861; Villumsen Sz Strauss 1987 : Harmon et al. 1987 ). optical catalogs (e.g. 



19871: 



Lahav et al 



Strauss et al 



1992 



19881: |Plionig||l988l). IRAS redshift surveys ( e.g., 



Webster et al. 



19971 ). and optical redshifts (|Lvnden-Bell et al 



Rowan- Robinson et al. 



1989; Hudso 



198a 
Lahav 
Il990|: 
udsonj 



19931 ). More d e tails about inhomogeneities on large scales in different kinds of surveys can be 
found in iLahavl (j200d ). where an extensive discussion is presented. 



An isotropically distributed population of radio sources is predicted to exist at cosmological 
redshifts in an isotropically expanding universe. The extragalactic radio sou rce population is ob- 



served to have a median redshift of 2 ~ 1 (e.g., iLoan. Wall. &: Lahavl 119971 ). a much later epoch 



than the cosmic microwave background. Assuming the radio source distribution is intrinsically 
isotropic, our local motion with respect to the rest frame of this population would produce a dipole 
anisotropy in the sky distribution of radio sources. Large-scale radio survey catalogs covering a 
large fraction of the sky can in principle be used to detect this dipole anisotropy and subsequently 
constrain the magnitude and direction of our local motion with respect to this frame. 

In this paper, I first describe the dipole anisotropy in the sky distribution of radio sources 
that would be introduced by our motion with respect to an isotropic radio source population. I 
then describe the corresponding dipole signal that would be seen in the distribution of discrete 
radio survey source counts. I describe a method to detect this signal in radio survey catalogs, with 
particular attention to the statistical significance of a possible detection. Finally, I consider both 
existing and proposed future radio surveys and whether a velocity dipole signal would be detectable 
with statistical confidence using this method given the number of discrete source counts present 
in these surveys. In this analysis, I consider only the most optimistic case in which there are no 
survey calibration problems, complications from the masking of certain survey regions, or uneven 
survey coverage. 



2. The Effect of Motion on the Observed Radio Source Distribution 



Motion with respect to the rest frame of an isotropic distribution of radio sources introduces 
two effects that change the observed source number density on the sky as a function of sky position 
([Ellis fc Baldwinlll984l : ICondonlll988l ) . The first effect is a boosting of the fiux from sources located 
toward the direction of motion. This depends on the observer speed and the source spectral index. 
This Doppler effect changes the number of detectable sources above a given fiux threshold in a radio 
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survey. The second effect, relativistic beaming, changes the apparent soUd angle on the celestial 
sphere as seen by the moving observer. This in turn alters the observed source number density as 
a function of sky position. For radio sources that are isotropically distributed across the full 
47r sr of sky, the source number density is (Tq = N/Att, a constant. The density variation from the 
observer's motion will depend on the angle 6 measured from the direction of motion (defined as 
= 0) so that (t{6) = dN/dVi. Azimuthal symmetry in the situation removes any dependence on 
the second angle (p. 

By combining the two effects described above, an expression fo r the source number d ensity 



of radio sources seen by the moving observer has been presented bv Ellis &: BaldwinI (| 19841 ). For 
small observer speeds (/3 = u/c ^ 1) in which only first-order terms in (5 are retained, the source 
number density function of 6 is shown to be: 



a(6') = {l + [2 + (7-l)(l + a)]^cos0}ao (1) 

which is a dipole anisotropy. In this expression, 7 is the power-law index for the number of 
extragalactic radio sources per flux interval {dN/dS ~ 5'^'^), and a is the mean radio spectral index 
of the source population (defined here according to 5 ~ i/^"). Following this, we can write: 



a{e) = (1 + ^cos6l)cro. 



(2) 



where the amplitude A = [2 -|- (7 — 1)(1 -|- a)](5. This expression (with some variation in 
nomenclature convention) has been used in the analys is of radio source counts by other authors (e.g. 



Condon 



198S; 



Baleisis et al. 



199^; 



Blake &: Wal]|l2002l ). The mean values of the radio spectral index 
and power-law index for the extragalactic radio source population are estimated to be a ~ 0.75, 
and 7 ~ 2, respectively (e.g. Blake & Wall 2002), and I use these values in the analysis throughout 
the rest of the paper. 



3. Detecting the Dipole Anisotropy in Radio Survey Catalogs 

In this section, I describe a method for detecting the dipole anisotropy in radio survey catalogs, 
and I calculate the minimum number of discrete catalog sources required to detect the anisotropy 
with statistical significance. For the calculation, I assume complete sky coverage (47r sr) in which 
the survey has no masked, incomplete, or invalid regions. Perfect survey flux calibration is also 
assumed. These conditions, of course, do not apply in real surveys, and I briefly address (but do not 
solve) the complications posed by the presence of these effects, which can be quite significant. Since 
survey masking and flux thresholds will reduce the number of usable source counts in a survey, I 
do consider the effect that this reduction has on the dipole detectability in two of the surveys that 
I discuss. 
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3.1. A Method for Detecting the Dipole Anisotropy 

A dipole vector D can be measured for a distribution of discrete sources by computing an 
evenly weighted sum over all source counts, where each source contributes a three-dimensional 
unit vector f to the sum (e.g., Baleisis et al. 1998)11 The direction of f is determined by the position 
of the source on the sky as seen by the observer at the center of the celestial sphere. The sum is 
then: 



N 



(3) 



1 = 1 



In the case of an ideal survey with perfect flux calibration, complete sky coverage, and in the 
limit of large N ^ all components of the vector sum would cancel except for a component along the 
direction of motion = 0), corresponding to the i-axis. 

In this case, the sum can be written as the integral of the number density over the celestial 
sphere: 



=27r p6=iT 

D = z I / a{6) cos 9 sin 
=0 Je=o 



(4) 



The additional cos 9 term is present here since we are only concerned with the z-component of 
the integral. Replacing cr{9) in the expression with Equation 2 gives: 



D = z 



=2n 



ct>=0 



which, when evaluated, gives: 



(70 (1 + AcosO) cos 6* sin 



(5) 



Att , 
D = — Aanz 
3 



(6) 



Using do = N/Ati and A= [2 + (7 — 1)(1 + a)]/3, the magnitude of the measured dipole can be 
written as: 



1^1 =Z) = ^[2 + (7-l)(l + a)]/?iV 



(7) 



^There is abundant previous literature on the use of directional statistics and spherical data and coordinates (like 
th ose employed here ) as applied ac r oss a numb er of scientific discipline s. Extensive discussions can be found in books 



bv lBatscheQ (|l98ll ). iFisher et all (Il993h . and iMardia fc JupdI j 19991 ) in which applications of directional statistics 
in a diverse range of subjects (e.g., biology, meteorology, psychology, earth sciences) are presented. 
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This is the dipole amplitude D one would measure for an ideal survey of N radio sources with 
complete sky coverage if the observer were moving with speed v = I3c (and v <^ c). 

3.2. The Statistical Significance of a Dipole Detection 

For a finite number of source counts A^, there will be an uncertainty in the measured dipole 
D arising from shot noise, and this uncertainty determines the significance with which D can be 
measured. We can imagine the uncertainty as being represented by a three-dimensional spherical 
probability cloud in direction space which represents possible measured values of D. This cloud 
will have some characteristic radial size, and it will be centered on the tip of the actual measured 
dipole D. The cloud is spherical since there is an equal probability of the deviation occurring in 
any three-dimensional direction. This probability cloud can be used to estimate the uncertainty in 
the measured speed of the observer /3, and exclude the zero-motion case (/3 = 0) at some confidence 
level (if an intrinsically isotropic distribution of radio sources is assumed) . The direction of motion 
can also be determined, with statistical uncertainties, by considering the direction of D and the 
angular size that the probability cloud subtends as seen from the origin (by the observer) . For the 
test of whether the zero-motion case is excluded with statistical significance, it is functionally more 
useful to re-center the probability cloud at the origin instead of centering it at the tip of D. This 
is done for the discussion below. 

The characteristic size p of the probability cloud is determined by a random walk of N steps 
in three dimensions (corresponding to the N sources in the vector sum that produces D). If the 
magnitude of each step is 1 unit (corresponding to the size of the unit vector) , then we must consider 
the average magnitude of these steps as projected onto an axis. The actual step sizes along this 
axis will range from —1 to with a probability weighting that is determined by the solid angle 
area from which the given projected value comes. The average projected step size s on the axis is 
computed from this to be s = l/\/3- 

For large N, the probability distribution of the magnitude of the final radial displacement is 

determined by the three one-dimensional Gaussian distributions in Cartesian coordinates which 
represent the three directional degrees of freedom. The probability p{r) that the displacement 
occurs at a radial distance r from the starting point (the origin) is: 

P^'^ = -^(i^2^^^' ^M-rV2p') (8) 

This is obtained by multiplying the three one-dimensional Gaussian probability distributions 
and including an additional weighting factor 47rr^ to account for the increasing spherical volume 
element as r increases. 

The characteristic size p is determined by the random walk according to: 
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p = sVn 




(9) 



Using this, the probabihty p{r) can be re-written in terms of N as 




exp(-3rV2iV) 



(10) 



Figure [T] shows the probabihty distribution p(r) from Equation [10] plotted against a histogram 
of the magnitude of the displacement vectors D that were calculated from 10^ simulated surveys. 
Each simulated survey had a random isotropic distribution of 2.7 x 10^ sources (which is the number 
of sources selected for the analysis and discussion of the combined 87GB/PMN survey below; see 
also Table [T]) . One of the simulated surveys is shown in Figure [2] and is discussed in more detail 
below. The histogram in Figure [1] was produced from the set of displacement vectors shown in 
Figure [3] (this is discussed below), and it has been normalized in Figured] to have a total area equal 
to one for comparison with the probability curve. The match between the curve and the histogram 
supports the use of this probability distribution in the subsequent calculations and analysis. 

We can use Equation [10] to estimate the likelihood of randomly measuring a dipole with 
a magnitude greater than or equal to D in the zero-motion, isotropic case. We integrate the 
probability function p{r) from r = D to an upper limit of r = (which is the maximum possible 
dipole value; in the limit of large N, the value of this upper limit approaches infinity). The 
integration is purely radial since the angular dependence drops out for the spherical cloud. The 
result of the integration, p(r > D), which ranges between and 1, is the probability- weighted 
volume outside of the radius D. This represents the likelihood of randomly measuring a dipole of 
magnitude D or greater from an isotropic source distribution. Conversely, we can use the measured 
dipole magnitude D to exclude the zero-motion case (/3 = 0) at a confidence level determined by 
1 — p{r > D) (e.g., if 1 — p = 68%, then the zero-motion case would be excluded at the la level). 

We can also determine an uncertainty in the magnitude of the measured dipole D using similar 
reasoning. We surround the measured vector D with the spherical probability cloud centered at its 
tip. We pick the confidence level of interest (Icr, 2a, etc.; let us call it na) from the corresponding 
percentage level of exclusion. Using this, we solve for the corresponding radial distance r„ in 
the probability cloud for which the probability-weighted volume within r„ is the aforementioned 
percentage of the total probability- weighted volume. This is done by evaluating the integral for 
p{r) using a lower limit (r^) which will give the appropriate percentage level of exclusion. Once r„ 
has been determined in this way, this defines an na error sphere centered on the tip of D. Since 



■^This result is not new, and the same expression can be derived differently by computing the probability 
distribution with three degrees of freedom, X3(3^^/-^)- Introducing a change in variables with a Jacobian and 
employing a substitution eventually leads to a result that is identical to Equation 1101 



-7- 



the uncertainty cloud is spherical, the projection of the cloud onto the D axis gives the magnitude 
r„, which is the na uncertainty in D. The measured dipole magnitude is then D zizrn (at the na 
confidence level). Figure H] shows detection confidence levels of dipole detections for calculations of 
D and r„ from several idealized surveys with different numbers of source counts N and different 
assumed values for the dipole velocity f3. 

This error sphere can also be used to determine an uncertainty in the direction of D. The 
projection of the error sphere onto the celestial sphere as seen by the observer at the center defines 
a circularly symmetric region on the sky centered on the dipole direction D. Sky directions lying 
outside the circular region are statistically excluded at the na level. Let the dipole direction D 
correspond to 6 = again, and let the na angular uncertainty in the dipole direction be 56n (this 
is the angular radius of the projected circular region). 56n is determined by the angle between D 
and a vector which starts at the origin and is tangent to the error sphere of radius r„. This is 
calculated by: 



for r„ < D. No statistically significant direction constraint can be made for the case r„ > D. 
This uncertainty angle 66n is determined by the number of sources A'' in the survey and the measured 
value of D (which in turn depends on N and the observer's velocity /3). Figure [5] shows calculations 
of 66n for an assumed dipole velocity of 370 km as a function of the number of survey source 
counts N. 

3.3. Incomplete Survey Sky Coverage and Imperfect Flux Calibration 

For real radio surveys which do not have complete sky coverage, an artificial dipole would 
be measured when the sum of the unit vectors f is taken over the observed sources. This 
dipole will be biased toward directions which have more complete survey coverage. To deal with 
this incompleteness when analyzing real radio survey data, one usually defines a mask of survey 
regions on the celestial sphere to exclude in the analysis. The vector sum is then taken over the 
distribution of sources in the unmasked, valid regions to obtain a measured dipole D. The center of 
the uncertainty probability cloud (described above) will be offset from the origin (corresponding to 
a magnitude and direction) by an amount determined by the geometry of the unmasked region. In 
principle, one can simply subtract this offset vector from all other vectors in the analysis (such as 
the measured D) so that everything (including the distribution of points in the probability cloud) 
is re-centered to the origin. One would then proceed with the analysis as described above. 

However, even without any problems introduced by the flux calibration, the shapes of the 
regions to be masked in a real survey are generally complicated. For example, sources in the 
Galactic plane and in certain declination bands may not be appropriate to include in the analysis 



69n = arcsin(r„/£') 
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(e.g., see the d ata m askin g used in the analysis of radio surveys by iLoan. Wall. &: Lahavl (|l997l ). 
Baleisis et aD (|l998l ). and Isiake k Walll (12003); see also Figure [2]). Also, the finite number of 
sources introduces a shot noise uncertainty in the offset vector. In the absence of calibration issues, 
the most straightforward course of action would be to simulate the survey by randomly assigning 
positions to N sources drawn from an isotropic distribution within the unmasked region of the sky 
and compute the sum of source unit vectors for many of these simulations. Figure [2] shows an Aitoff 
projection in equatorial coordinates of a single simulated survey of 2.7 x 10^ sources distributed 
isotropically over the unmasked regions of the sky (cf. Figure 1 of Baleisis et al. (jl998l ) and their 
masked regions, which were used as a guide here). The set of resulting vector displacements from 
the survey simulations will produce a probability cloud in which the mean position defines the 
vector to be subtracted. Figure [3] shows the vector displacements produced from 10^ simulated 
isotropic surveys like the one shown in Figure [2j In this plot, the mean vector has already been 
subtracted so that the cloud is centered at the origin, and the resulting displacements D have 
been projected onto a two-dimensional plane. This distribution of displacements in Figure [3] is the 
probability cloud that was used to produce the histogram shown in Figure [H 

Source clustering at local distances can also affect the analysis of radio surveys (e.g., the NVSS 
Survey; see discussions by Boughn & Crittenden 2002, Blake and Wall 2002, and also below). In 
this paper I make the assumption that the intrinsic distribution of sources is homogeneous and that 
there is no significant clustering on large scales beyond the local mass distribution. I assume that 
local sources would be excised prior to analysis. 

Also, real radio surveys do not have perfect flux calibration, and this can bias the source count 
distribution in certain sky regions, thereby producing an artificial dipole anisotropy. The effects of 
survey calibration errors in combination with the masking of survey regions can be quite severe 
and d ealing with these complex effects has been much discussed in th e literature (e.g. 



and 



(I2OO2I) for the ana l ysis o f the cosmic microwave background, and [Hamilton &: TegmarkI (120041 ) 



Hivonetal 



(York et al 



Swanson et al.l (j2008l ) for large-scale galaxy catalogs, such as the Sloan Digital Sky Survey 



2000)). In the analysis presented here, I ignore all such calibration effects and consider 



only the raw source counts N and whether a dipole would be detectable in the idealized case where 
the surveys and instruments are perfect. 



4. Prospects for a Dipole Detection in Existing and Proposed Future Radio Surveys 

Large-scale radio surveys currently exist in which the method outlined here for detecting the 
velocity dipole can be applied. However, the statistical significance of such a detection depends on 
whether these surveys have a sufficient number of sources, even if the flux calibration were perfect. 
Next-generation radio instruments that are in various stages of planning and development will also 
be used to conduct large-scale radio surveys in the future. Below I consider existing and possible 
future radio surveys in the order of increasing number of source counts, and I discuss the possibility 
of detecting the dipole in each of the surveys in turn, assuming the most optimistic scenario (i.e.. 
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perfect flux calibration). 



4.1. The 87 Green Bank and Parkes-MIT-NRAO Surveys 



The 87 Green Bank (H7GB) dOregorv fc CondoJl99lh and Parkes-MIT-NRAO (PMN) surveys 



(|Griffith et alJll994l : IWright et al 



1994 



19961') were conducte d at 4 .85 GHz and covered the nor thern 



and southern celestial sky, respectively. Baleisis et al. (jl998l ) and Loan. Wall. &: Lahav ( 1997 ) give 
concise descriptions and summaries of these surveys in the context of large-scale structure and the 
velocity dipole effect, and they refer t o a number of previo us attempts to detect the velocity dipole 
(see also the references in Section 1). iBaleisis et al.l (j 19981 ) have conducted an analysis of whether 
the 87GB and PMN surveys could be used for detecting the dipole using a method similar to the 
one described here (in their case, they considered the measured dipole magnitude relative to the 
shot noise term but did not do a more extensive analysis of possible detection significance) . 



In their analvsis. IBaleisis et alJ (119981 ^ excised certain incomplete or unreliable survey regions 
identified by Loan. Wall. Sz Lahav ( 199?! ). including sky regions within 10° of the Galactic plane, to 
eliminate Galactic sources, and they imposed a minimum flux density cutoff of 50 mJy for sources 
to be considered. This latter selection reduced declination-dependent number density variations. 
After performing a similar excision (but retaining only sources between 50 and 100 mJy), I find that 
the selected sample consists of ~ 2.7 x 10"^ sources covering ~ 70% of the celestial sphere (see Table 
[1]). Figure [2] shows similar mask ing used in a s i mulat ed survey of 2.7 x 10^ isotropically distributed 
sources. As was also found bv IBaleisis et al.l (j 19981 ). I find that this is not an adequate number 
of source counts for a statistically significant dipole detection given the expected local velocity of 
V ~ 370 km s~^, even if these surveys had perfect flux calibration. As seen in Figure SI an observer 
velocity that is an order of magnitude greater than this would be required for even a 3a detection. 
Figure [5] shows that for the expected velocity, no direction constrai nt could be made w ith this 
sample at any significance level. The shot noise analysis conducted by IBaleisis et al.l (|l998l ) showed 
that at least ~ 4 x 10^ galaxies over the sky would be needed to detect the velocity dipole at the 
same level as the shot noise (i.e., a la detection). Our results are consistent with theirs and suggest 
a comparable number of sources would be needed for a la detection (see Figure H] and Table [H 
particularly the NVSS select survey that has ~ 3.1 x 10^ sources). This result does not account 
for the dipole introduced from large-scale structure, which complicates the analysis further. Thus, 
the 87GB and PMN surveys are not adequate for detecting a velocity dipole anisotropy using this 
method, regardless of whether slightly different flux density selection criteria are used in the data 
selection. 
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4.2. The NRAO VLA Sky Survey 



The NVSS is a large-scale radio survey that was conducted with the Very Large Array at 1.4 
GHz ( Condon et alJ EoQS). The total number of sources in the NVSS is ~ 1.8 x 10^, and the survey 
is estimated to be 99% complete down to an integrated flux density of 3.5 mJy. The survey covers 
82% of the celestial sphere, corresponding to a declination range 6 > —40°, and the majority of the 
sources in the survey are believed to be at cosmological distances. 



Blake fc Walll (|2002l ) have searched the NVSS for a dipole anisotropy, and we follow their work 



to determine which survey regions and flux ranges could be considered reli able for the dipole sea rch 
described here. Regions within 15° of the Galactic plane were masked bv lBlake &: Walll (|2002l ) to 
eliminate Galactic sources in their dipole analysis. They also eliminated the "clustering dipole" 
(from the Local Supercluster), which has its own dipole component, thereby ensuring that only 
cosmological sources were used in the analysis. This was done by eliminating radio sources within 
30" of nearby galaxies k nown from several c atalogs. Although the claimed 99% completeness level 
of the NVSS is 3.5 mJy (ICondon et al.lll998l ). there are significant declination effects evident at this 
flux le vel in which the number density fluctuates by a few percent from the mean. iBlake fc Wall 
(|2002l ) use a minimum flux cutoff of 15 mJy for their analysis (see their Table 1), which reduces 
the variations to less than 1%. However, as described above, these flux variations in combination 
with survey masking may still introduce complications which can significantly contaminate a dipole 
search. 



After eliminating unreliable survey regions, iBlake &: Walll (|2002l ) retained a selected sample of 
~ 3.1 X 10^ sources above 15 mJy, representing ~ 20% of the initial sample (see Tabled]). This is 
not adequate for a statistically significant dipole detection with the method presented here, even 
if flux calibration were perfect in the survey. For an expected dipole velocity of 370 km s~^, the 
best possible dipole detection would be slightly less than la (Figure Hj). The dipole direction would 
also be unconstrained (Figure [5|). Given these issues, I conclude that the NVSS cannot be used for 
detecting the velocity dipole using this method. Including all of the sources in the entire NVSS 
catalog {N ~ 1.8 x 10^) would yield a marginal detection and constraint at best (~ 3a), but flux 
calibration problems with the weakest sources in the survey prevent this from being feasible. 



4.3. Next-generation Surveys with LOFAR and the SKA 



Future large-scale radio surveys will be conducted with next-generation radio telescope facilities 
that are currently in various stages of plann ing and deve lopment. Among these new instruments are 
the Long Wavelength Array (LWAfl (e.g., iTavloii 120061 ). the Murchison Widefleld Array (mwaJI 



^http: / /Iwa.unm.edu 

^http://www. haystack. mit.edu/mwa 
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(e.g., [Morales et al.l 120061 ) . and the Low Frequency Array (lofarA all of which will observe the 
sky at low radio frequencies (a few hundred MHz or less). Farther into the future, an even more 
advanced radio facility, the Square Kilometer Array (SKAjll, is being considered for development. 
Although the details of the surveys to be conducted with these instruments are not yet firmly 
established, I focus on two possible surveys that have been outlined for two of these instruments: 
LOFAR and the SKA. 

LOFAR is an advanced radio telescope array that is expected to operate at low radio frequen- 
cies (30 to 240 MHz) and will have t housands of antenna elements distrib uted over hundreds of 
kilometers (e.g., IStappers et al.l 120071 : iFalcke et al.l 120071 : iFender et al.l 120081 ). A survey using LO- 
FAR to search for gravitational lenses has been outlined by Jacksoi|§. In a possible LOFAR survey, 
half of the celestial sky (27r sr) would be covered down to a limiting sensitivity of ~ 0.7 mJy at 151 
MHz. The total number of sources expected to be detected in such a survey would be ~ 3.5 x 10^ 
(Table d]) . Other recent estimates of source counts from proposed LOFAR surveys at a variety of 
wavelengths suggest a range of ~ 10^ to ~ 10^ detectable sources in 2tt sr of sky (see, for example. 



the presentation by iRottgeringj 120071 ) . For the sake of simplicity, we use the survey parameters 
outlined by Jackson for our analysis in which the expected number of source counts falls near the 
middle of this range. 

The SKA is a planned next-generation radio telescope facility which will have vastly increased 
sensitivity for large-scale radio surveys. Simulations of a large-scale 1.4 GHz radio survey with 
the SKA suggest a lini iting flux density of ~ 0.1 /uJy and a yield of ~ 10^ or more sources per 
sr ( Hopkins et al. 1999 . [2OOQ ) . This corresponds to ~ 10^" sources for the entire celestial sphere 
(Tabled]). 

Even with large gaps in sky coverage, the sheer number of sources in a survey conducted with 
either LOFAR or the SKA would far exceed the requirements for a statistically significant dipole 
detection with this method. If flux calibration problems can be sufficiently minimized or properly 
corrected, and if local source contamination can be removed, it should be easy to detect the dipole 
with large statistical significance if the dipole velocity is near the expected value of ~ 370 km s~^ 
(see Figure H]). The constraint on the dipole direction should also be very precise (Figure [5]), easily 
constrained to within a few degrees or less in both cases. 



^http: / /www. lofar.org 

^http: / /www. skatelescope.org 

*N. Jackson, LOFAR Memorandum Series #4 



- 12 - 



5. Conclusions 

I have described a method for detecting a velocity dipole anisotropy in large-scale radio surveys 
and considered the feasibility of detecting this dipole with statistical significance in existing and 
proposed future large-scale surveys. This analysis does not account in any way for the severe 
complications that arise from imperfect flux calibration and masking effects in surveys, and therefore 
corresponds to the most optimistic detection case possible. Neither the combined 87GB/PMN 
survey nor the NVSS has a sufficient number of sources to detect the velocity dipole anisotropy with 
statistical significance using this method, even if no calibration or sample bias effects were present. 
However, proposed large-scale radio surveys using next-generation radio science instruments (e.g., 
LOFAR and the SKA) are more promising: surveys with these instruments should easily have 
enough source counts for a statistically significant dipole detection and direction constraint if flux 
calibration problems and contamination from local sources can be sufficiently reduced or eliminated. 

I thank Elizabeth Praton for assistance with some of the derivations, Nathan Keim for con- 
tributions to the simulation work, and Steve Boughn and the referee John Ralston for helpful 
comments and insights which have improved this work. 

REFERENCES 

Baleisis, A., Lahav, O., Loan, A. J., & Wall, J. V. 1998, MNRAS, 297, 545 
Batschelet, E. 1981, Circular Statistics in Biology (London: Academic Press) 
Bennett, C. L., et al. 2003, ApJS, 148, 1 
Blake, C. & Wall, J. 2002, Nature, 416, 150 

Boughn, S. P., & Crittenden, R. G. 2002, Physical Review Letters, 88, 021302 
Boughn, S. P., Crittenden, R. G., & Koehrsen, G. P. 2002, ApJ, 580, 672 
Condon, J. J. 1988, Galactic and Extragalactic Radio Astronomy, 641 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, R. A., Taylor, G. B., k Broderick, 
J. J. 1998, AJ, 115, 1693 

Ellis, G. F. R. & Baldwin, J. E. 1984, MNRAS, 206, 377 

Falcke, H. D., et al. 2007, Highlights of Astronomy, 14, 386 

Fender, R. et al. 2008, in Bursts, Pulses, and Flickering: Wide-field Monitoring of the Dynamic 
Radio Sky, Proceedings of Science (arXiv:0805:4349) 



-13- 



Fisher, N. I., Lewis, T., &: Embleton, B. J. J. 1993, Statistical Analysis of Spherical Data (Cam- 
bridge: Cambridge Univ. Press) 

Gregory, P. C. & Condon, J. J. 1991, ApJS, 75, 1011 

Griffith, M. R., Wright, A. E., Burke, B. F., & Ekers, R. D. 1994, ApJS, 90, 179 
Hamilton, A. J. S., & Tegmark, M. 2004, MNRAS, 349, 115 
Harmon, R. T., Lahav, O., & Meurs, E. J. A. 1987, MNRAS, 228, 5P 

Hivon, E., Gorski, K. M., Netterfield, C. B., Crih, B. P., Prunet, S., & Hansen, F. 2002, ApJ, 567, 
2 

Hopkins, A., Ekers, R., Jackson, C., Cram, L., Green, A., Manchester, R., Staveley-Smith, L., & 
Norris, R. 1999, Publications of the Astronomical Society of Australia, 16, 152 

Hopkins, A., Windhorst, R., Cram, L., k. Ekers, R. 2000, Experimental Astronomy, 10, 419 

Hudson, M. J. 1993, MNRAS, 265, 72 

Lahav, O. 1987, MNRAS, 225, 213 

Lahav, O. 2000, Cosmic Flows Workshop, 201, 377 

Lahav, O., Rowan-Robinson, M., & Lynden-Bell, D. 1988, MNRAS, 234, 677 

Lineweaver, C. H., Tenorio, L., Smoot, G. F., Keegstra, P., Banday, A. J., & Lubin, P. 1996, ApJ, 
470, 38 

Loan, A. J., Wah, J. V., & Lahav, O. 1997, MNRAS, 286, 994 
Lynden-Bell, D., Lahav, O., & Burstein, D. 1989, MNRAS, 241, 325 
Mardia, K. V. & Jupp, P. E. 1999, Directional Statistics (Chichester: Wiley) 
Meiksin, A., & Davis, M. 1986, AJ, 91, 191 

Morales, M. F., Bowman, J. D., Cappallo, R., Hewitt, J. N., & Lonsdale, C. J. 2006, New Astronomy 
Review, 50, 173 

Plionis, M. 1988, MNRAS, 234, 401 

Rottgering, H. 2007, in Astrophysics in the LOFAR Era, Emmen, Netherlands 
Rowan-Robinson, M., et al. 1990, MNRAS, 247, 1 
Smoot, G. F., et al. 1992, ApJ, 396, LI 



-14- 



Stappers, B. W., van Leeuwen, A. G. J., Kramer, M., Stinebring, D., Hessels, J. 2007, in Proceedings 
of the 363. WE-Heraeus Seminar on Neutron Stars and Pulsars, 100 

Strauss, M. A., Yahil, A., Davis, M., Huchra, J. P., & Fisher, K. 1992, ApJ, 397, 395 

Swanson, M. E. C, Tegmark, M., Hamilton, A. J. S., & Hill, J. C. 2008, MNRAS, 387, 1391 

Taylor, G.B. 2006, in Long Wavelength Astrophysics, 26th Meeting of the lAU, Prague, Czech 
Republic, 12, 17 

Villumsen, J. V., & Strauss, M. A. 1987, ApJ, 322, 37 

Webster, M., Lahav, O., & Fisher, K. 1997, MNRAS, 287, 425 

Wright, A. E., Griffith, M. R., Burke, B. F., & Ekers, R. D. 1994, ApJS, 91, 111 

Wright, A. E., Griffith, M. R., Hunt, A. J., Troup, E., Burke, B. F., & Ekers, R. D. 1996, ApJS, 
103, 145 

Yahil, A., Walker, D., & Rowan-Robinson, M. 1986, ApJ, 301, LI 
York, D. G., et al. 2000, AJ, 120, 1579 



This preprint was prepared with the AAS IM^^X macros v5.0. 



- 15 - 



0.008 



n I I I I r 



1 I I I I r 



0.006 - 



0.004 



0.002 



0.000 








100 



200 



300 



400 



500 



Fig. 1. — Histogram of the magnitude of the displacement vectors D calculated from 10^ survey 
simulations, each of which had 2.7 x 10^ sources isotropically distributed across the sky (see also 
Figure [3l from which this histogram was produced) . Also plotted is the probability distribution 
p{r) from Equation [10] as a function of radial displacement r, with = 2.7 x 10^ used in the 
probability expression. The histogram has been normalized in the plot to have a total area equal 
to one for comparison with p(r). The match between the curve and the histogram supports the use 
of this probability distribution in our analysis. 
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Fig. 2. — Aitoff projection plot in equatorial coordinates of 2.7 x 10"^ simulated radio sources 
distributed isotropically across the unmasked regions of t he sk y. T he masked regi o ns sh own here 
are very close to the ones used by iLoan. Wall. &: Lahavl (jl997l ) and iBaleisis et al.l (jl998l ) in their 
analysis of the combined 87GB/PMN survey (see also the discussion of these surveys in the text 
and Tabled]). The 2.7 x 10^ sources in this simulated survey correspond to the number of sources 
in the 87GB/PMN su rvey in the uri n iasked regions which have flux densities between 50 and 100 
mJy. See Figure 1 of IBaleisis et al.l (119981 ) for a qualitative comparison of this simulated source 
distribution with the actual source distribution. 
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Fig. 3. — Two dimensional projection plot of the set of displacement vectors D in three dimensions 
that were produced from 10^ simulations of isotropic surveys, each of which had 2.7 x 10^ sources. 
One of these simulated surveys is shown in Figure [2j The distribution shown here was produced 
after the mean of the distribution was subtracted from each displacement vector to correct for the 
vector offset introduced from survey masking effects. This three-dimensional distribution was used 
to produce the histogram of radial displacements shown in Figure [TJ 
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Fig. 4. — Dipole detection confidence level vs. dipole velocity for several existing and proposed 
future large-scale radio surveys. The solid curves represent the various surveys with different 
numbers of sources and are labeled (see the text and Tabled] for survey descriptions). In this plot, 
each survey is assumed to have perfect flux calibration and no contamination from local sources, 
corresponding to the most optimistic detection case possible. The expected dipole velocity of 370 
km s~^ is indicated by the dashed vertical line. For this velocity, a minimum of 2.0 x 10^, 3.1 x 10^, 
and 4.5 x 10^ survey sources would be required for 3a, 4c7, and 5a dipole detections, respectively. 
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Fig. 5. — Uncertainty in the measured dipole direction (69) vs. number of survey source counts 
(N) for several existing and proposed future large-scale radio surveys (see the text and Table 
[T]). In this plot, a dipole velocity of 370 km s^^ is assumed as is perfect flux calibration and no 
contamination from local sources (corresponding to the most optimistic detection case possible). 
The dotted vertical lines represent the various surveys with different numbers of sources and are 
labeled. From left to right, the solid curves represent la, 3a, 5a, and 8a constraints on the dipole 
direction uncertainty that are possible with this detection method. 
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Table 1. Existing and Proposed Future Large-Scale Radio Surveys. 



Radio Survey 


Observing Frequency 
(GHz) 


Flux Density Limit 
(mJy) 




87GB/PMN selecf^ 


4.85 


50 


2.7 X 10^ 


NVSS select^ 


1.4 


15 


3.1 X 10^ 


NVSS alF 


1.4 


3.5 


1.8 X 10*^ 


LOFAR'i 


0.15^ 


0.7 


3.5 X 10*^ 


SKA^ 


l.# 


10-4 


-lOio 



Note. — "Select" refers to a selected sample of sources from the survey deemed 
appropriate for the analysis. 

'^Includes only sources with a flux density between 50 and 100 mJy from the 
combined 87GB and PMN surveys. Several regions, including regions within 10 



degrees of the Gal actic plane, were also excluded (jLoan. Wall. &: Lahav 



1997; 



Baleisis et al.lll998l ). 



''Includes only sources with a flux density above 15 mJy from the NVSS. 
Several region s, including regions within 15 degrees of the Galactic plane, were 
also excluded (iBlake fc Walll l2002l ^ . 

'^Includes all sources in the NVSS with a flux density above 3.5 mJy. 
'^Includes all sources in a possible LOFAR survey outlined by Jackson. 
^Includes all sources in a possible SKA survey. 
^Possible observing frequency. 



